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ABSTRAC T 

Solutions  to  the  Dinchiet  and  Neumann  problems  for  the  region 
exterior  to  the  intersection  of  two  iegions  whose  individual  electro¬ 
static  Green's  functions  arc  known  are  developed.  The  method  is 
applied  specifically  to  obtain  solutions  for  the  extenor  of  a  solid 
finite  circular  cone  with  a  spherical  cap.  The  solutions  to  the  vector 
and  scalar  wave  equations  for  long  wavelength  can  be  expressed  in 
terms  of  these  Dinchiet  and  Neumann  solutions.  This  will  be  the 
subject  of  a  forthcoming  report. 
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SCOPE  OF  REPORT  AND  SUMMARY  OF  RESULTS 
Our  ultimate  purpose  is  to  obtain  solutions  to  the  scalar  'rave  equation 

v 2  0  +  k2  0  =  0 

over  the  exterior  of  a  compact  set  G  in  3 -space,  subject  to  vanishing  nor¬ 
mal  derivathes  on  the  boundary  of  Q ,  and  in  the  presence  of  a  source  ^  — 
i.e..  essentiall>  the  Green's  solution.  We  intend  to  get  explicit  solutions, 
moreover,  in  terms  of  classical  functions,  and  indeed  in  terms  of  the 
Green  s  functions  for  regions  simpler  than  In  this  paper  we  study  only 
the  case  k  =  0  (potential  problem),  and  we  will  subsequently  use  these  results 
to  obtain  the  so-called  Rayleigh  solutions  for  k  ^  0. 

This  program  is  accomplished  under  the  following  two  limitations: 

(1)  The  region  f2  is  the  intersection  of  two  regions  flj,  and 
for  each  of  which  the  exterior  electrostatic  Green's  function  is  known. 

(2)  The  region  Q  is  moreover  orally  symmetric  as  is  the 
function  C.  and  these  axes  of  symmetry  coincide. 

Using  the  present  method  it  does  not  appear  possible  to  lighten  these 
restrictions  cssentialh  (it  will  be  seen  under  (1)  that  we  could  have  an  arbi¬ 
trary  number  of  regions  instead  of  just  two).  These  two  restrictions. 
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which  correspond  roughly  to  the  two  major  parts  of  this  work,  are  imposed 
only  3s  they  are  needed  to  "make  the  method  work",  since  a  greater  generality 
at  the  intermediate  steps  adds  nothing  to  the  complexity  of  the  tieutment  and 
is  conceivably  of  independent  mterest. 

In  Part  I,  in  which  we  assume  only  the  restrictions  under  (1),  we 
solve  the  general  Dinchlet  problem  for  a  region  C  which  is  the  inter¬ 
section  of  two  regions  and  02  —  1.  a.  f  fi2  I°r  which  indi¬ 
vidually  the  electrostatic  Green's  functions  Gj  and  G2  are  known.  The 
method  is  probabilistic  in  nature,  using  Brownian  motion  theory,  ft  is  not 
indispensable  to  use  this  theory  to  ackeive  our  results,  but  it  was  the  way 
it  v*  as  discovered  and  it  has  the  advantage  of  being  relatively  simple,  and 
in  addition  gives  simple  bounds  on  the  iterative  method  it  leads  to. 

In  Part  II  we  assume  conditions  (1)  and  l2)  and  using  the  result  of 

Part  I  we  show  how  to  solve  Uie  Neumann  problem  (vanishing  normal 

v 

derivatives  on  the  boundary  of  ft)  with  a  suitable  singularity  (the  axially 
integrated  Green's  solution).  This  gives,  as  a  special  case,  ot  coarse, 
the  dipole  source  at  infinity  directed  along  the  axis  of  symmetry  analogous  to 
the  "radiation  condition”  at  infinity).  The  method  employed  is  due  essentially 
to  Bassett  and  is  very  old,  but  <«*ems  not  to  have  been  noticed  by  later  writers . 
the  present  modification,  in  particular,  seems  new. 
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I 

PROBABILISTIC  METHODS  AiiD  Til*  DUUCBLET  PROBLEM 

K  I  Notation  and  Terminology 

'fbe  points  of  3-space  E*  *re  denoted  by  vectors  p,  q,  \,  \ ,  , . . 
which  will  be  co-ordmati2ed  when  necessary.  We  consider  only  scalar- 
valued  functions  f,  g,  U,  V, . . .  The  volume  differential  at  o  is  denoted 
b\  d  v(p). 

3 

Wc  consider  certain  regions  £2,  £21#  £2^  ...  in  E  which,  if 

bounded,  w ill  be  compact.  The  complement  of  £2  is  denoted  by  £2  and 

the  boundary  of  U  bv  B  (£2).  We  denote  by  -  a  derivative  in  the 

3np 

direction  of  the  outward  drawn  normal  to  B(£2)  at  the  point  p  e  B(£2) 

The  restrictions  placed  on  £2  here  and  m  later  developments  will  be  clear 
from  the  context.  For  p  c  B(£2)  v  e  denote  by  d  a  (p)  the  surface  area 
differential. 

B\  the  electrostatic  Green’s  function  associated  with  £2  we  mean 

the  function  G(\,  >),  x,  y  c  £2  wh.ch  as  a  function  of  j  is  harmonic  in  fi 

except  at  \  =  a,  which  vanishes  on  B(£2),  and  is  such  that  G(x,  \)  ^ 

1 

-  ,  \  -*  If  B(£2)  is  considered  as  a  (possibly  unbounded) 

4r  »x  -  \ » 

grounded  conductor,  G(\,  \)  is  tne  potential  mduced  at  \  by  a  unit  charge 
at  \. 


4 


THE  UNIVERSITY  OF  MICHIGAN 

287 1 -5 -T 

1.2  Probability  Considerations 

We  let  \.(t),  0  <  t  <  co  be  the  three  dimensional  Wiener-E instem 
stochastic  process,  i.  e.  x(t)  =  £xt  (t),  x2  (t),  (t);  where  (t)  arc 

three  independent  one -dimensional  processes  with  x(0)  =  0  Let  a  compact 
Q  be  given  and  a  furction  f(p),  p  £  Bfft),  be  defined  over  B(f2).  Hie 
follow  m«  probabilistic  fac»  is  tb^  cornerstone  of  the  method  employed  here. 

Let  x0  €  Cl,  and  let  T  -  T(x0>  be  the  time  at  which  x(t)  +  xQ 
first  mtersects  Cl,  if  it  ever  does.  If  x(t)  +  \Q  never  intersect*  Cl, 
set  T  =  oo.  Thus 

Tfx^  =^supo  |t  |  x(r)  A  Q  . 

From  the  continuity  of  x(t)  we  conclude  x(T)  +  xQ  c  B (Cl)  and  we 
define  a  random  variable  =  V(xQ)  b\ 

\  f  (x(T)  +  xj)  ,  T  <  oo 

i 

i 

|  0  ,  T  =  oo 

Then  it  is  true  that  U(xQ)  -  E  (y(xJ)  ,  where  E  is  the  expecta¬ 
tion,  is  the  solution  to  the  exterior  Dirichlet  problem.  That  is,  U(x)  is 
haraiomc  for  x  c  Cl,  assumes  il;e  value  f(x)  for  \  c  B(C2),  and  U(\) 

vanishes  like  - — ,,  for  }xl  -»  oo. 

I\1  '  * 
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These  facts  are  relatively  well  known  —  cf  Ll,  2],,  j3j  —  and  may 
be  deduced  from  the  more  elementary  fact  that  if  Gq(x,  y )  is  the  electro¬ 
static  Green's  function  for  Q,  then  d  a  (p)  is  the  (infini- 

3np 

tesimal)  probability  that  x(t)  ■*-  Xq  first  enters  Q  via  the  surface  element 
d  o  (p)  at  p.  We  then  obtain  by  integration 


av(x0); 


J 

B(fi) 


f  (p) 


9np 


d  o  (p) , 


which  is  the  classical  solution. 


1 . 3  The  Dinchlet  Problem  for  Two  Intersecting  Regions 

Let  G  be  the  intersection  of  two  regions  and  C2  ~  Cj  Pfy 

and  suppose  that  the  electrostatic  Green's  functions  and  G2  are  known 

for  and  Q2  -  Exploiting  the  probabilistic  interpretation  of  the  Dinchlet 

problem  of  Section  1.2,  "e  show  ho  »  to  solve  the  Dirichiet  problem  for  C . 

The  method  will  .  Iso  yield  the  electrostatic  Green's  function  for  fi (which  we 

need  in  Parts  II  and  HI)  by  the  simple  expedient  of  solving  the  Dinchlet 

problem  for  those  boundary  values  induced  on  B(f2)  by  the  potential 

- — -  .  For  convenience  we  assume  at  least  one  of  the  two  regions 

4ff  I  \ 

and  Cl2  1S  bounoed,  so  that  the  same  is  true  of  f2  =  0  ^  ♦  7T»e 

method  is  applicable  if  H  is  the  intersection  of  Q,  ,  flj .  Qn  where 

the  individual  Creen’s  functions  are  known,  but  the  calculational  complexity 
increases  rapidly. 


6 


THE  UNIVERSITY  Or  MICHIGAN 
2871 -5-T 


Let  B{Oj.)  =5:  Us-/  an  1  B  £l2)  -  S*  U  Sy  v.here  Sj  =  , 

Sj’  -  B^jO  fk  and  similarlj  for  S2  ,  $2’ ;  S2  =  B(fl2)  0 
S21  =  B(ST22 )  P  fij.  To  make  the  problem  non-tnvial  we  assume  these  four 
components  non -void. 

In  connection  with  this  notation  we  refer  to  Figure  1  below  in  which 
o.  is  an  infinite  cone  and  fk  if  a  sphere,  an  example  treated  in  detail  in 
part  III. 
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Here  the  rachus  of  the  sphere  is  r  and  the  semi  vertex  angle  ol  the 
-1 

cone  is  s  -  cos  !uq. 

Let  f(p)  be  the  desired  boundary  value  defined  over  B(S2)  =  SjL  Sj 
and  consider  a  point  pc  ily  ifk.  Let  us  denote  by  Ejfp),  pc  SJ  the  required 
expectation,  which  is  the  potential  sought,  where  the  Wiener  process  starts 
at  a  point  pc  Si  and  similarly  Eo  (p),  pc  SJ  for  a  starting  position  pe  SJ  . 

For  the  particle  starting  at  pe  SJ,  it  may  intersect  Sj  before  intersecting 
SJ ,  or  SJ  before  Sx  ,  or  may  never  intersect  either.  If  we  denote  by  Gv  (p,  q) 
and  G2  (p, q)  the  respective  electrostatic  Green’s  functions  for  and  the 
sum  of  the  expectations  for  these  three  cases  is 

E,  (pi  =  F;  (p)  +  |  (p— "  E,  (x)aa(x),  (2) 

Si 

pc  S’  , 

2 

where 


Fj  (p)  = 


I 


3G,  (p,x) 
3n 

x 


f  (x)  d  a  (x), 

pr  Si  .. 


(3) 


In  an  exactly  similar  way  we  find  the  corresponding  equations  for  a 
point  pc  SJ 
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E,  (p) 


where 


<p)  +  . 


Si 


5G2  (p,  x) 
dn 

x 


t,  ix)  d  a  (x) 


pe  Si 


(4) 


F2  (p)  =  [  3  3^<1J,X)  Kxldo (x)  (5) 

I  X 

S2  p€  S} 

Equations  (2)  and  (4)  constitute  a  pair  of  litegral  equations  for  Ej 
and  E2  which  we  propose  to  solve.  We  remark  that  if  we  have  do«c  so,  and 
Fx  and  F2  are  defined  for  an  arbitrary  point  pc  ffjO isl2  by  (3)  and  (5;,  then 
the  right  hand  sides  of  either  (2)  or  (4)  gives  the  potential  slight,  and  they 
must  agree  xn  this  region.  If,  however,  as  in  Figure  1,  the  region  n2  is 
bounded  and  is  not  then  for  pc  Ujf  Vfiji  only  (4),  and  not  (2),  gives  the 
correct  result.  Similarly,  for  pc  Hz  »  (2)  only  gives  the  correct,  re¬ 
sult.  Indeed  it  is  easily  seen  that  the  above  enumeration  of  cases  is  valid 
for  both  (2)  and  (4)  only  if  pc  £2|Pi 
Let  us  define 

Hi  (p,q)  =  9.9. 1  -P*  qc  SJ  (G) 

Dn 

q 

H-.  (p,  q)  =  3^',p’q>  .  q  c  Si  n) 

a"q 
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*+  (P)  =  Fi  (p)  +  j  H,  (p.x,1  F,  (>)  do  (\) 

'  1  i  w 

j 

si  pc  a 

(p)  =  F,  (p)  +  ^  IU  <p,x>  F,  (x)doW  (g) 

./ 

a  p£3 

K,(p,q)»  II, 'P,  .\l  H2  (x,  q)  d  a  (\) 

J  (10) 

s;  p,  qc  a 

K,  (P.q)  --  J  H2(p,x)H1(x,q)do(x)  (n) 

S5  P,  qc  SJ 

"hcrc  F>  aaa  are  **“*«  in  (3)  and  (5).  Then  (2)  and  (4)  yield  the  follow, „K 
Fredholm  equations  for  E,  and  E2  , 


E,  (p, )  (p)  +  j  K,  tp,  x)  E,  (x)  d  a  (x). 

Si  s 


<P)  =3<P>  +  j  K,  (p,  x)  (x)  d  a  (x) 

s‘  pc  s; 

These  arc  two  ordinaty  Fredholm  equaUons  which  determined  E, 
and  E,  uniquely  and,  as  mentioned  above,  Uiey  completely  determine  the 


’  >  i-S. 
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i.4  1  he  Iterative  Solutions  for  Ei  and  E? 

In  these  integral  equations  the  classical  iteration  procedure  can  be 
carried  out  prolitably  (in  contrast  to  the  double -lajer  integral  equation  for 
which  the  iteration  procedure  diverges). 

(n)  th 

Considering  El  for  example,  we  let  Kj  vp,  q)  be  the  n  iterate 
of  the  kernel  Kx 
(1) 

Ki  (p»  q)  =  Ki  (p»  q) 

(n-M)  l  (n) 

Ki  (p.  q)  =  \  Ki  (p,  x)  Ki  (*>  q)  do  (x) 

J 

S2’ 

Then  \e  have  the  classical  Neumann  expansion 


\  i  i  % 

Ex  (p)  =  Fi  (p)  +  y  \  Kj  (p,  x)  Fx  (x)  do  (x) 

n  *  i  J 

V 


and  this  series  converges  exponentially  fast. 

To  verify  this  last  remark,  we  notice  that  (p,  q)  do  (q)  «s 
the  probability  that  the  Brownian  motion  particle  has  made  n  double 


il 


y  v  a  o  v 


m  n  i  j  ▼  t:  i  3  j  v 

'  •  '  ■  A  .  iv. 
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transfers  between  83'  and  Sj’,  terminating  at  an  area  element  da  (q) 
on  Sj’  w  ithout  ever  having  entered  Sj  p  S2  =  If  Hj,  say,  is 

bounded  it  is  clear,  at  least  intuitively,  that  the  event  of  n  such  double 
transfers,  for  large  n,  has  a  probability  0(an)  for  some  0  <■  a  1. 

More  precisely,  put 


a  =  sup 
x  c  Sz' 


K(.\,  q) 

’  V 


then  a  <  1  and  it  follows  easily 


Ki  (x,  y )  da  (y)  £  an  ,  n  -  1,  2,  . . . 


and  also  if  we  put 


b  =  max  J  Fi  (x)  | 

xe  .V 


we  shall  have 


E,  (p)  =  F,(x)  +  \  n 


Ki  (p,  q)  F:  (q)  •*  a  (q)  J 


* 

f- 


& 


> 


/ 

J 
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where 


I J  !  £  ^  ba1'1 

i  =  n  +  1 


b  a 


giving  ?  s,  Tuple  bound  on  the  eriur  made  by  truncating  the  series  at  n 
tc  rms. 

It  is  easily  seen  that  if  we  know  G2  (p,  q)  for  a  region  ft2  which 
,f  differs”  but  little  from  ft,  then  a  is  very  small,  and  as  a  matter  of  fact, 
pursuing  this  remark  it  is  possible  to  get  an  exact  perturbation  fonnula  for 
a  region  of  the  form  ft  =  ft2  U  where  T€  is  a  region  whose  capacity 
goes  to  zero  with  e. 

In  the  practical  calculations  we  have  encountered  thus  far  all  the 
iterates  can  be  explicitly  calculated  in  terms  of  matrix  calculations,  and 
the  resolvent  kernel  for  these  integral  equations  is  expressible  as  the 
inverse  of  an  explicit  matrix,  cf  Part  III. 
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SOLUTION  OF  THE  NEUMANN  PROBLEM  IN  TERMS 
OF  THE  DIRICH  LET  PROBLEM 


2  1  Assumptions  and  Notations 

From  this  point  on  we  assume,  as  stated  in  the  introduction,  that 
H  =  ft.  H  f2*  is  axially  symmetric,  or  is  a  body  of  rotation  The  sets  ^ 
and  ^2  need  not  be,  however.  We  also  adopt,  from  this  point  on,  the 
spherical  co-ordinate  system  (p,  0,  0)  to  describe  the  vectors  p,q, . . .  . 

In  Jits  notation  p  is  the  distance  from  the  origin  p  =  |  p|  0  is  the  co- 
latitudp  i'  £  0  ^  r  and  9  is  the  longitude  0  i  9  ±  2  ir .  We  adopt  the 
transformation/!  =cos  0,  so  that  1  >  -1,  and  consider,,  when  convenient, 

the  system  (p,/i,0). 

We  choose  the  axis  of  symmetry  of  S2  to  be  the  polar  axis  0=0.  Given 
a  region  12  we  say  that  a  function  U  (x),  x  t  Q,  is  a  solution  of  the  Neumann 
problem  if  U  (x)  is  harmonic,  x  e  12,  and  3/ 3n^  U  (x)  =  0,  x  e  B(f2).  In 
order  to  avoid  yielding  the  trivial  result  U  =  const,  the  condition  of  harmon- 
lcity  in  f2  must  be  relaxed  and  U  must  possess  certain  singularities  either 

3 

:n  the  finite  part  of  F  or  ‘at  infinity”.  The  precise  nature  of  these  is 
discussed  later 
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2. 2  Harmonic  Functions  and  Stream  Functions 
2  2 

We  denote  oy  ^7“  and  E  the  operators 

3  .  2  pi  .  5  2.  di  \  1  d^ & 

1  dp  3P  3p  V  *V  1-p2  ae2 


2  2 

2  2  a  e  ..2.  a 

e  g  =  p  — ?-  +  (1  -u  )  — 

dp" 


(15) 


We  apply  E“  only  to  axially  symmetric  functions,  so  the  term  in  6  does  not 
appear.  We  say  that  a  pair  of  functions,  axially  symmetric,  are  conjugate, 
if  they  satisfy  the  Beltrami-Stokes  equations.  That  is,  0  and  h  are  conjugate 
if  and  only  *f 


2  30  0  h 


(16) 


2  90 

P  9p 


9  h 

9p 


(17) 


It  follows  easily  that  if  0  and  h  are  conjugate,  0  must  satisfy 
0=0  and  h  must  satisfy  E^h  =  0.  Moreover,  if  0  and  h  are  conjugate 
they  are  orthogona1  -  V  0  *  yh  =0.  The  pair  of  equations  (16;  and  (17) 

correspond  to  the  ordinary  Cauchy-Riemann  equations,  and  given  any  v 
2 

satisfying  y  j  =0  there  corresponds,  oy  means  of  them,  a  function  g 

2  2 
satisfying  E  g  =  0,  and  conversely;  given  g  satisfying  E  g-0  there  ex.sts 


■  >  .  i-h. 
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2 

a  0  satsifying  Vj  0=0.  The  correspondence  is  unique  up  to  an  addluve 
constant  which  maj  be  determined  by  prescribing  the  value  at  one  point,  or 
at  infinity  We  give  the  exact  reciprocity  relations  later 

The  goal,  which  we  described  in  the  introduction,  of  solving  only 
Dinchlet  boundary  value  problems  could  be  achieved  if.  given  n,  we  could 

o  _ 

solve  the  problem  E“h  =  0,  p  €  o  h(a)  =  const,  a  e  B(ft),  for  then  the  func¬ 
tion  0  conjugate  to  h  would  be  harmonic  with  the  Neumann  boundary  conditions 
It  turns  out  to  be  possible  to  do  this  by  solving  an  accessory  Dinchlet  problem 


2.  o  ihe  Method  ot  Bassett 

Basrett  notirprf  that  the  following  relation  is  true,  for  anv  function 
u  =  u(p,  p)  with  requisite  derivatives. 


cos©  E 


v/l-p2  u  (p,  p)^  =  p  y/l-ii2  V  L  ^cos$  u(p,  .  (18) 


This  correspondence  between  E  and  V7 »,  winch  is  easily  proved  by 
substituting  in  (14)  and  (15),  is  the  link  by  which  wTe  connect  the  Dinchlet 
and  Neumann  problems.  Bassett's  results  are  credited  to  him  by 
Hobson  in  (¥] .  A  related  procedure,  due  to  Weinstein  and  Payne,  is 
sketched  briefly  In  Part  V. 

From  (18)  it  is  clear  that  if  we  want  to  solve  the  problem  E2g  =  0  with 
boundary  values  b(p),  d  €  B(f2)  we  can  do  so  by  solving  ^2V  =0  with 
boundary  values  c(p>  =  -p=  b(p),  and  in  fact  the  function  g  wll  then  be  given  by 

*1“M^ 
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g 


'/  1-u 


COS0 


(19) 


If  we  now  fmd  a  function  0  conjugate  to  g,  0  will  be  harmonic  and  its 
equipotcntial  surfaces  will  be  orthogonal  to  the  level  surfaces  of  g 

In  particular,  suppose  v(x)  is  harmonic  in  ft,  except  for  prescribed 
singularities.  In  order  to  solve  the  Neumann  problem  we  consider  the 
function  V  conjugate  to  v  and  choose  b(p)  =  v(p),  p  e  B(Q).  Then  the 
function  v  -  g,  with  g  as  above  has  zero  boundary-  values  and  its  conjugate 
function  v  -  t  is  harmonic  and  has  vanishing  normal  denvates  on 
and  thus  famishes  the  solution  with  the  prescribed  singularities  This 
procedure  is  in  fact  the  way  the  further  development  proceeds. 

We  remark  that  the  restriction  to  axially  symmetric  bodies  and 
boundary  conditions  is  indispensable  for  the  use  of  (18j,  and  within  this 
framework  cannot  be  extended 


2  4  Solution  of  Reciprocity  Equations 

Wc  consider  the  equations  (16)  and  (:7),  showing  how  to  express  0  m 
terms  of  h  and  vice  versa. 

Suppose  0  iu  an  axially  symmetric  harmonic  function  defined  or.  fi. 
Then  there  exist  measures  and  v2  over  the  non-negative  reals  such  that 
for  sufficiently  small  p 
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0  ip.  u)  -  j  p*P  (m)  di/x  (X)  *20) 

and  for'  p  sufficiently  large 
o 

0  (p,  p)  =  ^  px<fj)  di,'!  w  •  (21) 

These  a.e  the  "normal  solutions"  and  generally  r,  and  r2  will  be  discrete 
so  that  the  above  integrals  will  be  sums  -  that  is,  uv  for  example,  will 
assign  a  weight  at  the  value  X^  ^  0,  n  =0, 1, .  .  . 

When  0  is  expressed  m  this  fashion  its  conjugate  0  =  h  has  a  particularly 
simple  form  We  obtain,  namely,  by  solving  (16)  and  (17),  that  for  p  in  the 
range  where  (20)  is  valid 

a  =  h  (p,  p)  =  -  p  /l"-pZ  ^  <p)  dt  j  (X)  (22) 

and  for  p  in  the  range  of  (21) 

0  =  h  (p,  *x)  =  p  /l-p2  ^  pJJW.GI  .  (23) 

These  formulas  determine  h  except  for  an  additive  constant.  We  omit 
the  derivation  winch  is  straif  htforward. 

Conversely,  any  axially  symmetric  stream  function  h  car.  be  written  m 
the  form  of  (22)  oi  (23)  w:ih  suitable  v 1  and  v2,  and  the  harmonic  function 
which  is  conjugate  to  it  is  given  by  (20),  (21)  with  these  same  i  x ,  i2 
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2.  5  Solution  of  the  Neumann  Problem 

Let  G(p.  p)  be  an  axially  symmetric  function  harmonic  in  fi  except 
for  prescribed  singularities.  Using  the  preceding  results  we  are  no*  in 
a  position  to  solve  the  Neumann  proolem  for  £2,  that  is  finding  a  function 
U  (p,  p)  harmomc  in  £2  with  vanishing  normal  derivatives  on  Q  and  with 
the  same  Singularities  as  G. 

We  write 

\J( p)  =  G(p)  -  V(p)  (24) 

so  that  V(p)  is  harmonic  everywhere  in  12,  and  we  require  it  to  vanish  icr 
p  — *  cu,  thus  requiring  (J  and  G  have  the  same  limiting  behavior. 

Write  G  as 

G  =  P^(p)  dr  (X)  ,,  P  4.  PQ  ,  (25) 

or.  supposing  the  singularities  of  G  arc  exterior  to  the  sphere  P  =  PQ  - 
The  function  conjugate  to  G  is  then,  by  the  preceding  section, 

C  pX 

G=-p  '/l-p2  ^  T^TT  dp(X) 

and  *e  write  for  a  point  p  e  B(£2) 

b(p)  =  -  ^  -^y  (p)  dr  (X)  cos*  .  (26) 
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Let  V\(p)  cgs0  he  the  solution  to  the  Diriclaet  problem  for  C  with 
boundary  values  b{p)  -  it  is  simple  to  prove  that  the  solution  is  of  the  form 
V,(p  pjcosd.  Then  p  V.(p,  p)  is  a  stream  function  and  Vj  can  thus 

be  written  tn  the  form  of  (22)  or  (23) 


Vj(p,  P) 


(p)dr,00 


V,(p,  p> 


P(xl)  (p)  dt2  (W 


(27) 


We  now  set  V  as  the  function  conjugate  to  p 
have  from  the  reciprocity  formulas  of  Section  2. 4 

V  =  -  ^  (X+l)  p*  P^(p>  dTj  (» 

v  *  5  p>(")dr’<» 


Vt ,  so  that  we 


(28) 


for  small  and  large  p  respectively.  The  requirement  that  V-*0,  p-*co 
then  shows  the  additive  undetermined  constant  must  be,  in  fact,  zero. 

Since  G  -  p  Vl  is  r  stream  function  vanishing  on  its 

con.ugate  has  normal  derivatives  zero  on  B(f2),  as  required. 

To  recapitulate,  the  normal  solution  to  the  Neumann  problem,  given 
n  and  G,  is 
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(the  bar  representing  a  conjugate)  where  \\  cc so  is  the  solution  to  the 
Dinchiet  problem  for  Q  with  boundary  values  b  given  by  (26),  y  being 
determined  by  (25).  The  function  V  =  p  Vl -p-  V!  is  given  by  .'28)  where 
Tj  and  X2  are  determined  by  (27) 


2  (.  The  Green’s  Solution 


We  now  specialize  the  solution  of  the  preceding  section  to  the  case 
where  the  singularity  is  that  generated  by  "charged  ring".  This  will  gi\e 
tne  Green’s  solution,  or  fundamental  solution. 

Thus  in  equation  (24)  we  take  for  the  function  G 


G<po  “o 


i 


end  this  can  be  expanded 


oo 

po'VP>p)  '  L.  ~^T  WPnW  •  p‘po 

^  Po 

CD 

P  (^  )  P  (|l)  p  >  p 

rH-1  no  n  o 


and  this  is  the  Green’s  function  for  the  operator  defined  by  (14)  with 
respect  to  the  differential  dp  dp.  That  is 
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n00  f  1 

V.  \  \  G(p  ,f!  ;p,|j)  S(p,p)dpdp  =  c(p  p  ; 

1  Jq  *>-I  °  °  o  o 


)  (31) 


for  continuous  functions  f 

Bj  (22)  and  (23'  fne  function  conjugate  to  this  G  is 
co 

G(po.po ;  p.p)  •  -  p  \^7  JZ  T^i  p‘UM  Pn(p0> 

1  (n+i)p 


oo  n 

G(po  po;p.p)  =  P  V'l  -p2  XZ  “ii 


The  term  n  =  0  does  not  appear  since  p'  (p)  =  0. 

By  (26;  we  consider,  lor  (p,p)  e  B(  «)  and  f or  p  <  p  when 


(p.p)  e  B(  52),  the  boundary  value 


b(p,p.  0)  - 


<n+l)r 


—  p[\)  Pn(pD)  cose 


and  we  let  V,(p,p) cosh  be  the  solution  to  tm  ordinary  Dinchlct  problem 
with  these  boundary  values  on  B(fl).  The  function  V,(p.p)  can  then  be 
expressed  by  (27),  in  the  normal  form 


4 

Winn 


mm 


I  i  *  A  1  3  J  V  .  S 


-yr  v  a  o  * 
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v^p.v)  =  ^  '  anf.  “  P*1* 1 

0 

00 

v.(p.p)  ■  )  -yV  pfi 
-  PA«  An 


for  small  and  large  p  respectively,  the  constants  depending  or.  p^ ,,  . 

For  the  function  conjugate  to  V1  we  have  by  (28) 


oo  > 

V(p,p)  -  ^]an(Xn+l)p"  Px 


q>  . 

-r. 

— *  "  „Xn 


fr  V* 


and  consequently  the  required  solution  to  the  Neumann  problem  is 


U  P»p)  '  G(p  .  *P  >  p)  “  V(p,  p) 

v  o  o  o  o 


where  G  is  given  by  (32)  and  V  by  (36) 

To  recapitulate,  given  we  solve  the  Dii  ichiet  problem  with 

boundary  values  b (p,  p  ,  6)  given  by  (34)  The  solution  ca  i  lie  expressed 

as  ViCOs8  which  determines  die  constants  a  =a  (p  n  )  in  (33).  These 

n  n  o  o 

determine  V(p.u)  via  (36),  and  then  U  is  given  by  (371 
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III 

AN  APPLICATION  TO  A  SOLID  FINITE  CONE. 

3. 1  The  Problem  Considered 

Using  sphencal  coordinates  let  ft  be  the  intersection  of 

U,:  P  <  rD 

The  relevant  regions  and  surfaces  are  as  »n  Fig.  1. 

Iu  this  part  u-e  are  going  to  carry  thi  ough  some  of  the  calculations 
for  ft  that  we  developed  in  pai  ts  I  and  II.  The  calculations,  though  elementary, 
become  somewhat  involved,  but  fortunately  the  use  of  matrix  notation  permits 
a  significant  compression. 

We  solve  here  the  Neumann  problem  for  ft,  corresponding  to  a 
singularity  at  infinity  of  the  form  py,  p — >  oo.  This  is  done  in  reasonable 
detail  and  the  answer,  given  by  (45),  is  expressed  in  terms  of  the  quite  sim¬ 
ple  matrix given  by  (43). 

For  the  Dinchlet  problem  for  the  exterior  of  ft,  that  :s  the  solution 
«ith  boundary  values  one  on  B(ft),is  similar,  except  that  it  is  not  necessary 
to  find  the  conjugate  functions,  and  is  given  by  (4fi). 

3. 2  The  Neumann  Problem 

As  mentioned  in  the  preceding  section  we  do  not  use  the  singularity 
generated  by  (30),  but  only  its  limiting  form  when  p<j— *  oo  andpQ  - 1. 
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This  will  bo  the  analogue  of  the  '  radiation  condition’’  and  is  that  induced 
by  uii  indefinitely  large  charge  on  the  remote  axis  of  symmetry  -  1.  e. 
The  potenu?‘  ueliavcs  pa.  P — >oo  as  can  be  seen  by  considering 


lim  p  (G  -  —  ) 
o  p 

p->  00  o 


m  equation  (30) 


By  considenn;  (32)  this  induces  the  boundary  values 


l:m  p^b  (p,*i,9)  =-  4pPx(J  (p)  cos  9 

P-*oo 


=  y  p  Vi  -f*2  cos  e 


when  we  set  pq  -  1  in  (34)  and  take  the  limit.  Thus  in  (26)  we  have  just  one 
term  and  the  boundary  value  for  the  associated  Dinchlet  problem  is 


b  ( p,n ,  0)  =  o  V  1  -  p“  cos  9. 


3. 3  Solution  ol  the  Associated  Dinchlet  Problem 

Let  (m)  be  the  positive  zeros  of  (p0)  considered  as  a  function 
ot  \#  j  =  1, 2, . . . ,  m  -  0, 1, . . .  and  put 

_  ar(m,M  I 


3P(m,(n  )  J 

X  ( m )  =  - —  I 

Ox  |  x  -  a,(m). 
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Then  (he  Greens  fi-ncuon  (or  (2,  is 

to  oo  a.  (m) 

-  2  (2  -  tfm,  0 )  pp 


4 (P. i) =  v_  -1TL1 

S^frr  < 1  - "o  > 4 *  £J <m>  + 1  “j  (m)  p) <m) 


cos  m  (Op  -  6^) 


if  p^  >  pp  ,*  if  Pq  <  pp  we  interchange  p  and  q  in  (39). 


(39) 


G?  (p,q)  r  ~r— 
■4  r 


For  ((3  ,  the  sphere,  we  have 

oo  oo  2n+l 

,  fo 
(1 


r»-T) 


2n  +  1  2n+  1 


(40) 


p!m)  <V)  P(nra,  <u  )  cos  m  (6  -  0  ). 
n  p  n  q  P  q 

for  0  -<  n  .  We  mane  a  similar  interchange  of  p  and  q  if  p  §.  p  . 
P  q  r  P  q 

For  B  ( 12,)  and  B  (S2.)  we  have 


d  a,  (p)  =  p  V  1  -  p  2  dp  d9 


S  1  J, 2 

3np  p  V  P0 


9 


Y  pc  B(f2,) 


'  "  =  "o  J 


d  o2  (p)  =  r“  d  Gdp 


a 

dp 


P  =r 


j 

< 

j 


pc  B (12,) 
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In  calculating  Hx,  H2.  Q.  Kj,  K2of(6),  (7),  (8),  (9),  (10) 
and  (11)  we  need  only  retain  the  terms  corresponding  to  m  =  1  in  (39)  and 
(10)  because  o i  the  boundary  condition  (38).  cf.  also  the  remark  following 
Uo). 

We  let  a  ,  a  ,  /3  be  respectively  a  (1),  o  (1),  0  (1),  and  in  the 
)  )  J  3  3  J 

sequel  we  replace  the  cosines  by  complex  exponentials,  it  being  understood 
that  in  the  ensuing  calculations  we  are  to  take  real  parts  and  subsequent  to 
any  0  integrations  we  are  to  supply  a  factor  r. ,  since  an  easy  calculation 
shows 


This  notational  device  enables  us  to  write  Hi  and  h2,  as  well  as 

the  other  relevant  functions,  as  rather  compact  bilinear  forms. 

Thus  retaining  only  the  terns  rn  =  1  in  (39),  (40)  and  defining  the 

following  column  vectors,  with  p  -  (p,#i,  b), 

/a  (1)  \ 

U,  (p)  --(p  1  Pa  <p)  e  j 


v ' 1  - 


(pi 


f^L 

\  a:  +  1 


*m  m  “l  i  J  ▼  a  a  T  S  *  .  s 
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-  /  p"1  1  £_ 

x  J  o  J 


uc  hutv  for  Ilj,  denned  b>  (6) 


Hj  (p,q)  -  UJ  (p)  Vx  (q)  ,  pp  <  Pq 


=  UJ  (p)  V,  (q)  ,  np  >  p^ 


qe  B  («,), 

where  the  accent  denotes  the  transpose  and  where,  without  specifically 
indicating  it  real  parts  are  to  be  taken. 


Similarly  defining 


)-inZ±  r  j- V1’  fc)e-!fl 
\  4tt  o  }  ) 


we  obtain  from  (40)  and  (7) 


Ii2  (p,q)  =  Y*  (p)  7.  (q) 


qe  B \il2) 


The  follow  two  vectors  A,  33  and  three  matrices  M,  N*:,  N2,  whose 
components  are  constants,  are  needed.  In  the  designation  of  a  matrix  as 
c  ,  the  index  i  refers  to  the  row  index  and  k  the  column  index.  In  denoting 
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bv  0  a  column  vector,  d  "efers  to  the  component  of  it. 
J  J 


fv 


A  ~  J  V2  (\)  b(x)da1(«> 
Si 


.  2  a.  +  2 

1  -  li  r  i 
O  O 


2£  (aj  +  2) 


B-  j"  Z  {x)  b  (x)  da,  (\) 

s' 


8  (j  -  1  )  (j  +  2) 


\ir. 


pii,cu  >.\T7 


O  J  o 


«.j  v, 


r3 

IT"  <3%  ■  "o  '  2  !  •  J 1 


(x)  Y'  (x)  dd!  (x) 


Pk>fco> 


(a  l)r  J 


+  k+  1 


J  =1.2. . 

k  =  0, 1, . . . . 


i%'t  =  J  Z(x)  UJ  (x)  da2  (x) 

s; 


A  . 
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=  ( 


1  \ 
-V-n  * 


2]  +  1 
4 


ak  +  -'+l  „  (1),  . 

ro  P, 

(a^  -  j  )  (a^  +  j  +  1  ) 


1=0.1, . 

k  =  1 , 2, _ 


X,  -  J  Z  {> 


(x)  Lf2  (x)  d  o2  (x) 


_7,  2  2  j  +  1 

-(1  -m0  >  — — 


p0)  K)av 

j  o  k 


(ak'J)(ak+J  +  1,r0 


ak  -  J 


j=0,l . 

k  =  1,  . 

In  these  integrations  the  function  b  (p)  is  given  by  (38)  and  a  factor 
r.  has  been  supplied  in  lieu  of  the  «  integrations,  as  explained  above. 

In  terms  of  these  expressions  it  is  simple  to  calculate  the  quantities 
in  (3)  and  (5),  and  (8)  thru  (11).  They  are,  namely,  the  real  parts  of  the 
left  hand  sides  of  the  following 


F,  (p)  =  Ci  (p)  A  ,  p  e  (  )  O2 

F2  (p>  -  Y'  (p)  B  ,  p  >  r 
P  0 


(p)  -  U$  (p)  A  +  UJ  (p)  M  B, 
(?)  =  Y*  (p)  (B  +  X2  A) , 


pe  s; 

pe  Si 
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Ki  (p,  q)  -  UJ  (o'  M  Z  <q) ,  p,  q  e  S£ 


K2  (p,  q)  =  Y'  (pi  N,  V,  (q)  ,  p,  q  c  SJ 
In  terms  of  these  expressions  wv  can  easily  solve  (12)  and  (13). 
From  (13)  we  obtain  for  E2  namely, 


E2  (p)  =  Y‘  (p)  B  t  Y>  (p;  s,A 


+  Y’  (p)  K,  I  V,  (x>  E2  (x)  d  o,  lx) , 


Denote  by  X  the  vector 


r 

X  -  V, 

V 


(x)  E2  (x)  d  c,  (x). 


then  multiply  (42)  by  V,  (p)  and  integrate  over  SJ ,  obtaining 


X  =  MB  i  M K2  A  1  M  K,  X 


X  -  (1  -  M  X’i )  (  M  E  +  M  N2  A  ) 

Substituting  this  in  (42)  we  obtain 


E2<p)=Y*(p;  [fl  +  XjA+X^I-MX,)"1  M(B  +  K2A)  j 
=  Y'  { p)  j ^  I  +  Nj  (  I  -  M  Nj)  "  AmJ  (  B  +  X2  A) 


-  Y*  <p)  (I  -  X,  M)~  (  B  -  X,  A  ) 
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;  tii?  vector  -/V  bv 


j\-  ( X  )  =  (  I  -  i\;  M)  1  (  B  +  N2  A  ) 


PU>  (k) 


which,  though  initially  defined  only  for  y  =u  (  p  c  S!  )  is  now  seen  to  be 

o 

valid  everywhere  in  p  >  (p  e  il^)  by  continuati?  a  this  being  a  consequence 
of  the  remark  fol lo^ng  (5).  This  remark  is  also  verified  after  we  calculate 

Ei. 

To  calculate  Ex>  we  obtain  from  (12) 


Ei  (p)  =  Ui  (p)  A  +  L*1  (p)  M  B 

f 

+  b;  <t<>  M  J  z  lx)  E,  (xl  d  o  lx) 
Si 


This  is  solved  exactly  as  for  E,,  and  tve  obtain 


Ei  (pi  -  l’i  (p)  A  -  l’i  (p)  M  B 


T  "1  !p)  M  (I  -  Nj  M  )  (X2  A  +  N,  M  B) . 


On  using  {4)  we  obtain  for  the  potential  W(p)  ,  p  t  o. ,  the  cxprcssio 


J2 
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♦«  H  V.T 


W  (p)  =  Y1  (n)  H  •  Y*  (pi  j  X,  A  +  N;  M  B 


-  X,  Mil-  X,  M  )  *  (  X,  A  +  X,  M  B  >  j 

=  \  '  (p)  B  (i  -  K,  M)  <X,  A  +  X,  MB  >1 
=  Y’  (p)  (J  -  X,  Ml  ‘  1  (X,  A  +  B  ) 

--Y'  (p.'A 

in  agreement  with  {43} 

A  similar  argument  with  (2)  yields 

f 

V  (?)  -  F,  'p>  +  J  II,  (p,  x)  E2  (x)  da  (x) 

s; 

=  UJ  (p)  A  -  J  US  (p)  Y;  (x)  Y’  'x)  da  (x) 

L  <\.<pp 

*  U|  (p)  |  V,  (X)  Y>  (x)  d  a  (x)  Hi-NjmJ*1 

'\>PP  J 
th 

For  use  j-.  component  of  the  xector  m  brackets  v»e  have 

„«> ®  , ....  r 


(B  +  Nj  A)  . 


,0  Pk  V  V7  %,_<»*>  fl _  _J _  1 

k  -  1  ]  a  -  k  f 

P  ‘  j  L  I  >  J 
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237 1 


'  s  .  *,  - 1  - 


Pf1!  IpJ  <-  *>  ‘  k 

k  co 


p  (3  ■  k) 
J  J 


J  =1 

m  which  we  have  set  p  ~  (  p,p,  0). 

In  the  first  term  the  summation  on  i  can  be  performed  and  the 

result  is  simply  p^(p)  /p(l)  ,  .  .  This  follows  directly  «n  expanding 
K  /  k  ^  o 


p[^  (p)  as  a  Fouriei  series  in  the  complete  orthogonal  setjp  Oi)  rover 

la‘  J 

the  interval  i  >  n  ^V-Q  >  the  sum  of  the  senes  equaling  this  function  except 
perhaps  at  the  end  points  of  this  interval. 


Thus,  defining  the  matrix  Q  by 


Q  = 


(1)  a  -  k 
P.  V)  r  J 
k  o  o _ 

Wk) 


J  =  1.2 . 

k  -  0, 1, . . . 


we  obtain  for  p  t 


•»  (Pl  =  UJ  (pl  A  +  Y'  (p)  j  I  -  N,  M  1  (B  +  N2  A) 

-  I' (  (p!  Q  f  I  -  N,  M  (BA  N,  A  ). 

I  — 

These  two  expressions  for  W,  nameiy  this  latte i  and  (44)  must  agree 
over  their  common  domain  of  definition,  i.e  f) 

The  potentiil  sought  is  then 

0(1), 

P  <**} 

W(P)=  ^  X.  -rTT-  cos9,p>ro. 
where  the  A.  are  the  components  ol/Y  given  by  (43) . 
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The  stream  function  for  the  body  f.  of  Fig.  1  is  then,  by  (23), 


- 2  V  ^ 


and  the  solution  to  the  Neumann  problem  for  ft  Is  by  (24)  and  (28) 

P„  M  (45) 

X(p)=pp-\  nXn  —  j  . 

1  '  P 

To  recapitulate,  let  X^  be  the  0th  component  of  the  vector  V\  , 

given  by  (43)  ,  and  A  by  (41).  Then  the  solution  to  the  Neumann  problem  for 

ii,  which  behaves  as  p/i  for  p-jro  ,  is  given  by  (45>. 

vYc*  remarl,  that  in  the  matrix  Q,  if  u  is  such  that  a  "  k  for  some 

o  j 

’Ptegers  ;  a"d  k,  that  element  in  Q  becomes  indeterminate,  and  is  to  be 
replaced  by  its  limiting  value  when  k-}a.  ,  re  rarding  k  as  a  continuous 
variable.  This  limiting  value  is  readily  seen  to  be  unity. 

Also,  as  explained  in  section  1. 4,  the  inverse  in  (43)  can  be 
approximated  by 


-  1  *  “  7 

(1  -  N»  M)  =  (\\  M)  A  J 


where  the  bound  for  J  is  known,  as  ir»  section  i.  4. 
n 


♦ 

‘5 

_  ....  =  •  i  kL _ •' _ : _ ; _ 

- 

- 

* 

N: 
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"V 

- 

*  -i  Solution  of  the  Conductor  Potential  Problem 

We  sohe  here  the  Dirichlet  piobiem  forO,  in  which  the  function 

- 

b  <p*  =  1  instead  of  ( 58),  and  seek  the  solutions  which  vanish  as  p  *  for 

p — ■)  oo.  The  calculations  are  formally  nearly  identical  to  those  in  the 

preceding  section,  except  that  because  of  the  simple  constant  b(p) 

_ 

the  0  inteip-ations  can  be  avoided 

- 

Define 

a  =  l**1  positive  zero  of  P  (p  ) 

J  x  o 

“j  3v  | 

3  P  (i  ) 

"j  9*  |  x=a 

and  the  following  vectors  A,  B,  and  matrices  M,  N1#  N2 

r  JJT  1 

-  ■  ~  ■ 

A  =  ° 

(a  +  1) 

J  3 

„  -  2l*1rMl  <-1>PS^o> 

■  » . 

D  2  ‘o  j  ( j  +  1)  '  1  '  ° 

4  ro‘wo+1>-  0 

- 

P!>o> 

■>  /  <■  a,  -p  k  +  i 

P  (a  +  k  +  1)  r  J 

1  1  o 

i 

.it. 

- 

( 

/ 

/ 

•tt  -  '  ■- 

. 

*  ' 

/* 
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ak  ~  J  +  1 


,  ,  2  2  j  +  1  o 


Pfc  >a, 

jo  k 


_W!l - - 

(\  '  j)  <ak  +  J  U  ro”k  ' 


These  correspond  to  the  previous  designated  quantities  with  the 
same  letters,  and  we  have  for  the  solution  of  the  Dirichlet  problem 
P  (m) 

D(p)=  >  X  —  — 

/  v  n  n  +  1 

*  O 

where  A  is  the  n**1  component  of  given  by 


-A.=  (i-n1m)'1<B+1 
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IV 

THE  METHOD  OF  BASSETT 

4  i  Preliminary  Remarks 

The  success  of  the  present  method  rests  very  intimately  on  Bassett's 
identity  (18)  .  This  identity  enabled  us  to  revert  the  solution  of  the  Neumann 
problem  to  that  of  a  Dinchlet  problem  There  is  thus  a  certain  interest  in 
attempting  to  extend  and  simplify  it,  if  possible 

In  this  part  of  the  report  we  give  a  sketch  of  the  method  of  Weinstein  and 
Payne  which  is  related  to  Bassett’s  method,  and  show  that  Bassett's  identity 
is  essentially  the  only  one  of  its  kind 

4  2  The  Uniqueness  of  Bassett's  Method 

The  crux  of  the  method  of  Bassett  consists  in  finding  a  universal  boundary 
function,  given  in  (26),  such  that  the  solution  to  the  Dinchlet  problem  corres¬ 
ponding  to  it  could  be,  under  a  suitable  transformation,  changed  to  the  Neumann 
solution  by  considering  the  conjugate  function  to  the  transformed  function  For 
die  method  to  he  universally  applicable  to  every  ft  the  boundary  function  and 
transformation  employed  must  be  independent  of  ft 

From  this  last  remark  it  is  simple  to  deduce  th«l  ihe  boundary  function 
must  be  a  harmonic  function  We  only  have  to  consider  a  sequence  of  boundaries 
ftn  converging  to  ft  and  employ  Hamack's  theorem  Also  the  transformation 

38 
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must  be  l.near  in  the  sens*-  that  *f  D  is  the  solution  to  tae  Dinchlct  piwuiem 
the  transform  of  D  mast  be  a  L>  +  b  where  a  and  b  are  certain  point 
functions 


If  we  let  aD  +  b-F  we  obtain  D  =  a F  +  £  for  some  functions  a  and  jj, 
and  *f  we  want  F  =  const  for  a  point  on  B(ft)  for  every  ft,  then  a  =  B  is  the 
universal  bounding  function  sought  and  £  =  const,  and  we  may  as  well  assume 
p  =  0  We  now  consider  only  axially  symmetric  ft  -  if  uniqueness  can  be 

established  for  this  subset  of  all  ft  it  will  follow  over  the  la*  ger  class 

2  2 

The  condition  required  now  is  that  V  BF  =  0  implies  E  F  =  0  and 

conversely  From  the  fact  that  the  classes  of  harmonic  functions  ami  stream 

2  2 

functions  »orm  linear  manifolds  we  have  then  S7  BF  =  CE  F  for  every'  F,  and 
we  seek  functions  B  and  C  such  that  this  is  an  identity  in  F. 

The  calculation  becomes  simpler  in  cylindrical  co-ordinates,  and  (14) 
(15)  become  with  co-ordinates  (z,r,  0) 

V  3z2  3r2  r  3r  r2  3«2 


so  that 


E2  -A 


3r 


I  JL 

r  dr 


2  o  n  l 

v  — 

r  3r  2 


90 


ii  we  are  then  to  have  the  identity  =  CE  F,  we  obtain  for  F  axially 


symmetric 
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V  BF  =  F  V“  B  *  2  ly  T  V  B  +  By2  F 


ce2f  =  cV"F  - 

r  dr 


If  these  expressions  are  to  be  identical  in  F  then 


V2b  =  0 


uhich  verifies  the  previous  remark  that  the  bounding  function  is  harmonic 
The  equation  in  F  then  becomes 


i  VF  \7  B  =  -  —  J>E_ 


which  shows  B  cannot  depend  on  z  Since  in  that  case  V  F  V  B  =—  ~ 

3r  dr 

we  obtain  that  B  must  satisfy 

dB  _  _B_ 

dr  r 

Q 

and  B  7  "+lore  c  is  some  function  of  9  Finally  if  B  is  to  hr  harmonic 
we  must  have  C  'tO)  *  C'(6)  -  0  and  C  =  A  eoso  -  A  sins 

Thus  we  obtain  B  =-  (A  cos8  +  A'  sin8)  or  expressed  in  spherical 
co-ordi  nates 

d  _  A  cos9  +  A*  sinfl 


!or  some  constants  A,  A'  The  choice  used  in  (19)  was  A  =  1,  A’  =0 
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4  t  Tne  Method  ot  Weinstein-Payne 

Using  what  they  *erm  the  "method  of  generalized  electrostatics”, 
Weinstein  and  Pavne[7J  have  developed  a  procedure  similar  to  the  above  to 
‘solve  the  stream  function  equation  in  terms  of  t?»c  solution  of  a  Dinchltt 
problem  Though  their  method  is  not  rele\ ant  to  the  development  presented 
here  it  seems  worthwhile  to  point  out  the  similarity 

Their  method  depends  on  the  fact  there  exists  an  identity  similar  to 
Bassett's  *n  the  case  when  the  dimensionality  of  the  problem  is  stepped  up  by 

2  9 

two  Using  the  cylindrical  co-ordinates  again  we  let  V  and  E~  be  the 

n  n 

Laplacian  and  stream  opeiators  in  n  dimensions  with  axial  symmetry  Thus 
for  n  >  ? 


The  identity  which  corresponds  to  {18)  is  then 


and  the  reciprocity  relations  (lf>),  (17)  are 


30  _  I  3c  _  30 _ 1 _ 3c_ 

Oz  n-2  dr  '  3r  n-2  3z 
r  r 
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Choosing  n  =  3  it  is  thus  seen  that  by  solving  the  conductor  potential 
problem  {the  electrostatic  problem)  in  5  dimensions  we  can  find  the  stream 
function  in  3  dimensions  and  thus  the  Neumann  problem  in  3  dimensions  This 
method  has  the  advantage  that  the  Dirichlet  problem  considered  is  one  with 
constant  boundary  values,  but  it  has  the  serious  disacu a.ntage  of  requiting 
5  dimensional  potentials  In  our  work  it  has  seemed  preferable  to  use  only 
3  dimensional  potentials  with  non-constant  boundary  values 
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